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Abstract 

We study J/ip photoproduction via fragmentation at next-to-leading order in 
the QCD-improved parton model, using the nonrelativistic factorization formalism 
proposed by Bodwin, Braaten, and Lepage. We consider direct and resolved photo- 
production of prompt J/ip mesons and Xcj mesons radiatively decaying to J/ip + j, 
taking into account the formation of both color-singlet and color-octet cc states. 
Adopting the values of the long-distance color-octet matrix elements extracted from 
fits to prompt- J/^ data recently taken at the Fermilab Tevatron, we predict that 
measurements of J ftp photoproduction at DESY HERA should show a distinctive 
excess over the expectation based on the color-singlet model at small values of the 
inelasticity variable z. This is complementary to the expected enhancement at z <, 1 
due to the color-octet contribution to photon-gluon fusion. 
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Some time ago, experiments P] at the Fermilab pp collider Tevatron have revealed that 
the production rate of prompt charmonium at large transverse momentum (px) exceeds 
the most accurate theoretical predictions within the so-called color-singlet model (CSM) 
by more than one order of magnitude. This surprising observation may be interpreted 
by taking into account two recent developments that have revolutionized the theoretical 
description of heavy-quarkonium production in high-energy collisions. Firstly, it has been 
realized |§ that fragmentation production must dominate at large px, which implies that 
most charmonium in this kinematic domain is produced by the hadronization of individual 
high-pr partons. Secondly, a rigorous factorization formalism for quarkonium production 
based on nonrelativistic quantum chromodynamics (QCD) has been developed [|J, which 
allows for a systematical treatment of the formation of charmonium from color-singlet 
and color-octet cc pairs. This formalism comprises a separation of short-distance parts, 
which are amenable to perturbative QCD, from long-distance matrix elements, which 
must be computed through lattice simulations or extracted from experiment, and it takes 
into account the complete structure of the charmonium Fock space. An important new 
feature of this approach is the appearance of the so-called color-octet processes, in which 
a cc pair is produced at short distances in a color-octet state and subsequently evolves 
into physical (color- singlet) charmonium by the nonperturbative emission of soft gluons. 
In fact, the Tevatron data [|IJ on prompt charmonium production at large px can be 
successfully interpreted by including the appropriate color-octet processes and adjusting 
their long-distance matrix elements to fit these data ||. 

In order to convincingly establish the phenomenological significance of the color-octet 
mechanism, it is indispensable to identify the color-octet contributions in other kinds of 
high-energy experiments as well. On the theoretical side, a first step in this direction has 
recently been taken in Ref. || by investigating the signature of color-octet charmonium 
production in e + e~ annihilation on the Z-boson resonance. Although the experimental 
data from CERN LEP @ are in agreement with the predicted signature, they do not 
exclude the hypothesis that the observed prompt-charmonium signal is produced by color- 
singlet processes alone. The analysis of low-energy e + e~ data might be more instructive 

In this somewhat unsatisfactory situation, it would be very desirable if the experiments 
at the DESY ep collider HERA could probe the color-octet mechanism of charmonium 
production and independently measure the magnitude of the color-octet matrix elements. 
HERA is presently operated in such a way that E e = 27.5 GeV positrons collide head-on 
with E p = 820 GeV protons in the laboratory frame, so that approximately 300 GeV 
are available in the center-of-mass (CM) system. Charmonium production at HERA 
dominantly proceeds via photoproduction, where the incident positrons (or electrons) are 
scattered by small angles and effectively act as a source of quasi-real photons, which 
interact with the constituents of the incoming protons via hard scattering. The scattering 
can either be elastic, diffractive, or inelastic. In elastic scattering the protons stay intact, 
while in diffractive proton dissociation they are destroyed by a small momentum transfer 
from the photons. Both mechanisms may be interpreted by assuming that the protons 
interact with their environs by the exchange of so-called pomerons, so that an adequate 
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theoretical treatment lies beyond the scope of the ordinary QCD-improved parton model. 
For a dedicated study of the color-octet processes, it is thus necessary to focus attention 
on inelastic scattering. In terms of the inelasticity variable z = p p ■ pj/^/p p ■ p 7 , where 
p p , p 7 , and pj/^j are the proton, photon, and J /if} four-momenta, respectively, this may 
be achieved by discarding events with z values close to unity. As may be inferred from 
the relation z = 2E p rriT exp(—yi ah ) /W 2 , where is the J/ip rapidity in the laboratory 
frame (with y\ ah > in the proton flight direction), m T = {M 2 ,^ + p\) 1 / 2 is its transverse 
mass, and W is the 7p CM energy, inelastic events tend to have large y^ and large W . 

Recently, the HI || and ZEUS collaborations presented their 1994 data on the 
inelastic photoproduction of J/ip mesons, collected in the kinematical ranges 0.45 < z < 
0.9, 30 GeV< W < 150 GeV and 0.3 < z < 0.8, 60 GeV< W < 130 GeV, respectively. 
After dividing out the respective integrated photon flux factors, these data were compared 



10 1 with a next-to-leading-order (NLO) calculation [11] of 7 + p — * J/ip + X via 



7<7 fusion |T2|] within the CSM, where W was fixed at the representative value W = 
100 GeV. Although reasonable agreement was found || [IIJ, one should bear in mind that 



the theoretical prediction strongly depends on the chosen input values for the charm- 
quark mass m c , the asymptotic scale parameter A, the renormalization scale /x, and the 
factorization scale Mf, with the normalization uncertainty being as large as a factor of 
three [[□]] . One might hence conclude that there is still some room left for alternative J / ip 



production mechanisms. 

Very recently, the analysis of Ref. [II] was complemented JR] at leading order (LO) by 
including the 2 — > 2 color-octet processes 7+ g — > cc[ 8, n]+g and 7+g — > cc[8, n]+q, where 
n = 1 Sq, 3 Si, 3 Pj with J = 0,1,2. Adopting the color-octet matrix elements determined 
from fits to the Tevatron data |l[ on prompt J / ip production, the authors of Ref. [|13 



found that the channels 7 + g — > cc[8, n] + g with n = ^o, 3 Po, 3 P2, the cross sections of 
which are divergent in the limit z — > 1, due to g — > gg collinear splitting, dramatically 
increase the cross section of inelastic J/ip photoproduction at HERA in the upper z 
range, at z^0.8. On the other hand, the color-octet contribution rapidly falls off for z 
decreasing, and is negligibly small already for z^0.6 ||13|| . 

Being higher twist, jg fusion is only relevant if the produced J/ip mesons have trans- 
verse momentum p T <^Mj/^ w 3 GeV. In analogy to the situation at the Tevatron 0, at 
Pt 3> Mj/^, one expects fragmentation to be the dominant mechanism of inelastic J/ip 
photoproduction at HERA, i.e., single partons (mostly c, c quarks and gluons) that are 
produced with high p? by the hard 'jp scattering fragment into J / ip mesons. The purpose 
of this Letter is to present a NLO analysis of this mechanism in the framework of the 
QCD-improved parton model, and to assess the potential of HERA to probe the color- 
octet processes in the range of low to intermediate z values, thus providing an independent 
check of the Tevatron measurements. 

A first study [[L4]] of inelastic J/ip photoproduction via fragmentation at HERA only 
included the Bethe- Heftier process 7 + g — > c + c and the Compton process 7 + q — » g + q 
followed by c — » J/ip and g — » J/ip fragmentation at LO. These processes contribute to 
direct photoproduction, where the quasi-real photons couple directly to the quarks in- 
volved in the hard scattering. In a fraction of time, the quasi-real photons fluctuate into 
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bunches of quarks and gluons, which in turn interact with the proton constituents via 
hard scattering, while the photon remnants give rise to hadronic activity in the backward 
direction (resolved photoproduction). The longitudinal- momentum distributions of the 
partons inside the resolved photons are described by parton density functions (PDF's). 
At NLO, the direct- and resolved-photon contributions both strongly depend on the fac- 
torization scheme and scale, and must be combined in order to give a well-defined physical 
observable. From previous experience in connection with inclusive pion, kaon and 
D*± |L6], [T7| photoproduction at HERA, one may infer that, in the case of J/ip photo- 
production via fragmentation, the resolved-photon contribution, which was disregarded 
in Ref. fl4| , is likely to be much more important than the direct one. Furthermore, one 



expects a substantial enhancement due to the NLO corrections, which were neglected in 
Ref. WM- 



The factorization formalism developed in Ref. 0] implies that the charmonium frag- 
mentation functions have the general form 

D^ H (x,fi) = J2di-> n (x,v)(0\O H [n}\0), (1) 

n 

where di^ n (x, /x) gives the probability for the parton i to form a jet that includes a cc pair 
in the state labelled by n carrying the longitudinal-momentum fraction x, and (0\O H [n] |0) 
measures the probability for a pointlike cc pair in the state n to bind to form a physical 
charmonium state H. The coefficient di^ n (x, /io) at the initial scale /i = 2m c ~ Mj/^ 
involves only momenta of order m c , and can thus be calculated in nonrelativistic QCD as 
a perturbation expansion in the running coupling constant a s (/i ) @- The evolution of 
Di^H^x, Ho) up to higher fragmentation scales /i is ruled by the timelike Altarelli-Parisi 



equations, which may be conveniently solved at NLO in x space |T7 . 

We take into account prompt J/ip mesons as well as Xcj mesons (J = 0,1,2) that 
radiatively decay to J/ip + 7 with well-known branching fractions |L8| (non-prompt J/ ip 
mesons). The dominant J/ip (Xcj) Fock states are [1, 3 Si] and [8, 3 Si] ([ 1, 3 Pj] and [8, 3 Si]), 
respectively. The scaling of the respective matrix elements with the mass m c and the 
relative velocity v of the bound c and c quarks is indicated in Table [TJ. In the case of J j ip 
(Xcj), the color-singlet matrix element is related to the (derivative of the) nonrelativisitic 
radial wave at the origin and may thus be extracted from the measured leptonic (light 



hadronic) annihilation rate pi. The color-octet matrix elements have been determined 



from fits H to Tevatron data [T| . The input values JT9[ adopted in our numerical analysis 



are also listed in Table |1|. In the case of D g ^j/^, the v A suppression of (0|O J// ^[8, 3 «Si]|0) 
relative to (0|O J /^[!, 3 5i]|0) is numerically compensated by the fact that d^rgag^ has 
two powers of a s less than d g ^i-i Sl y Such a compensation does not occur for D c ^j/^ and 
D c _> XcJ , so that color-octet fragmentation is negligible in these cases. The color-singlet 
and color-octet contributions to D g ^ XcJ are theoretically intertwined, since d g ^[i, 3 Pj] has 
a logarithmic infrared singularity, which is absorbed into (0|O XcJ [ 8, 3 5*i] 



We work at NLO in the MS scheme with nt = 4 massless quark flavors, A^- — 
296 MeV [^tj, and \i = Mf = tut- We take c and c to be active partons inside the proton 
and resolved photon, for which we adopt the CTEQ4M |20| and GRV \21\ PDF's, respec- 



tively. As described in Ref. [16|], we adjust the factorization of the final-state collinear 
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singularities associated with c and c in such a way that it matches the massive calculation. 
We treat the quasi-real photon spectrum in the Weizsacker- Williams approximation with 
a maximum virtuality of Qmax = 4 GeV 2 as described in Ref. [16|. We compare our frag- 
mentation results with the LO prediction for direct J /if) photoproduction via 7$ fusion 
within the CSM 0, which we evaluate with the QCD-corrected value of (0\O J ^[ 1, 3 Si]|0) 
specified in Table [3]. In this way, we include the bulk of the NLO enhancement, the resid- 
ual dynamical QCD correction being of order 20% [[llj]. Direct Xcj photoproduction via 
7<? fusion is forbidden at LO in the CMS, and it is marginal in the color-octet channel |22| . 
Resolved J /if) photoproduction mediated by gg fusion, which receives both color-singlet 
and color-octet contributions, visibly contributes only at very low z and does not change 
the qualitative picture of jg fusion [p2fl . 

From Fig. we observe that the fragmentation mechanism vastly dominates inelastic 
J /if) photoproduction in the lower z range. For a minimum-p^ cut of 4 GeV (8 GeV), its 
contribution exceeds the one due to jg fusion for z < 0.4 (0.75), by factors of about 4 
and 200 (20 and 700) at z — 0.25 and 0.05, respectively. The bulk of the fragmentation 
contribution is induced by resolved photons, except for z close to unity, where fragmenta- 
tion is anyway insignificant. Since we wish to elaborate color-octet signatures at low and 
intermediate z, we regard jg fusion, which for z ^0.6 is well described by the CSM |i"3| , 
as a background. As is evident from Fig. |2], we may efficiently suppress this background 
by introducing a stringent maximum- z cut, of 0.4 say. In fact, if the maximum- z cut is 
lowered from 0.8 to 0.4, the cross section of 7^ fusion at px = 2 GeV (10 GeV) is reduced 
by a factor of 4.2 (2.8), while that of fragmentation is only insignificantly decreased, by 
11% (15%). At p T = 2 GeV (10 GeV), 2% and 30% (4% and 24%) of the fragmentation 
cross section integrated over 0.05 < z < 0.4 are due to the prompt color-singlet and 
non-prompt channels, respectively. Thus, this regime offers a unique laboratory to probe 
the prompt color-octet channel and thus to measure (0|OW[8, 3 5i]|0), so as to check the 
Tevatron result. Even for 0.05 < z < 0.8, the fragmentation cross section in Fig. ^| is in 
excess of the one due to jg fusion. As may be inferred from Fig. [j], this is subject to 
change if the minimum- z cut is increased. This is illustrated in Fig. |3|, where the ZEUS 
J /if) — > data jnj collected in the interval 0.3 < z < 0.8 are compared with the 



respective fragmentation and 7g-fusion predictions. We have divided the experimental 
data points by the estimated extrapolation factor 1.07 which was included in Ref. JII| to 
account for the unmeasured contribution from < z < 0.3. In fact, our combined analysis 
of fragmentation and 7^ fusion suggests that, at pt = 2 GeV (5 GeV), this factor should 
be as large as 5.2 (3.4). The corresponding value for 0.05 < z < 0.3 is 2.6 (2.6). The gap 
between the combined result for 0.3 < z < 0.8 and the ZEUS data is partly filled by the 
dynamical NLO corrections |Tl| and the color-octet contributions Jl3[ to 7g-fusion, which 
are not included in Fig. |[ 

In conclusion, the cross section of inelastic J / if) photoproduction in ep collisions at low 
z and large px is very sensitive to the color-octet matrix element (0|(9 J// ^[8, 3 Si] |0). We 
propose to accordingly extend previous measurements of this cross section at HERA, in 
order to obtain an independent, nontrivial check of the Tevatron color-octet charmonium 
puzzle. 



5 



Table 1: Values |T9| of the matrix elements used in the numerical analysis and their mass 
and velocity scaling. The relations between the J-dependent matrix elements follow from 
heavy-quark spin symmetry. 



(OpW[T^p 1.13 GeV 3 

(0|C> J ^[8, 3 Si]|0) 0.014 GeV 3 

(0\O x * J [l, 3 Pj]\0)/(2J+ 1) 0.0880 GeV 5 

(0|C9^[8, 3 S' 1 ]|0)/(2J r + 1) 0.0076 GeV 3 



[m 3 v 3 ] 
[m 3 v 7 ] 

[m 3 t> 5 1 
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Figure 1: Cross section da/dz of inelastic J/ip photoproduction at HERA, integrated 
over 60 GeV < W < 280 GeV and 4 GeV < p T < 12 GeV (upper curves) or 
8 GeV < pt < 12 GeV (lower curves). The NLO fragmentation contributions due to 
direct photons (dotted lines), resolved photons (dot-dashed lines), and their sum (solid 
lines) are compared with the LO 7g-fusion contribution (dashed lines). 
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Figure 2: Cross section dcr/dp T of inelastic J/ip photoproduction at HERA, integrated 
over 60 GeV < W < 280 GeV and 0.05 < z < 0.8 (upper curves) or 0.05 < z < 0.4 
(lower curves). The total NLO fragmentation contribution (solid lines) is compared with 
the LO 7g-fusion contribution (dashed lines). For comparison, also the prompt color- 
singlet (dotted lines) and the non-prompt (dot-dashed lines) fragmentation contributions 
are shown. 
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Figure 3: Cross section da/dp? of inelastic J /if) photoproduction at HERA, integrated 
over 60 GeV < W < 130 GeV and 0.3 < z < 0.8 (lower curves) and divided by the 
corresponding photon flux factor, 6.66 x 10 -2 . The total NLO fragmentation (solid lines) 
and LO 7^-fusion (dashed lines) contributions are compared with the ZEUS data ]IT|. 
For comparison, the theoretical results are also shown for 0.05 < z < 0.8 (upper curves). 
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